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H 要 : 近年 来 黄河 流域 人 与 水 资源 矛盾 愈 发 突出 ,尤其 是 地 下 水 储量 的 过 度 消耗 已 成 为 限制 该 区 


域 人 与 自然 和 谐 发 展 的 主要 矛盾 之 一 。 基 于 重力 反 演 与 气候 实验 卫星 数据 和 全 球 陆 盏 


数据 同化 


系数 数据 ,根据 流域 水 循环 与 水 量 平衡 原理 ,测算 2003 一 2016 年 黄河 流域 地 下 水 储量 的 变化 情况 ， 
在 探究 其 时 空 变化 特征 的 基础 上 ,识别 了 地 下 水 储量 下 降 明 显 区 域 的 人 口 暴 露 度 。 结 果 表 明 :(1) 


空间 分 布 上 ,黄河 流域 地 下 水 储量 西 多 东 少 ,由 


向 东 下 降 程 度 不 断 加 剧 , 且 下 降 区 域 呈 现 由 下 游 


向 中 上 游 扩散 的 态势 。 下 降 区 域 主要 集中 在 中 东部 地 区 ,变化 等 级 以 较 剧烈 减少 和 剧烈 减少 为 
主 ,(2) 时 间 变 化 上 ,2003 一 2016 年 黄河 流域 地 下 水 储量 整体 呈 下 降 趋势 ,年 均 降 幅 为 5.93 mma 
其 中 ,2004 一 2016 年 为 连续 下 降 ,2015 一 2016 年 的 下 降幅 度 最 大 ;另外 ,地 下 水 储量 时 间 变 化 存在 
季节 效应 ,地 下 水 储量 下 降幅 度 春季 最 大 冬季 次 之 ,秋季 再 次 之 、 夏 季 最 小 。(3) 地 下 水 暴露 风险 


状态 下 的 地 级 市 人 口 密度 在 空间 上 呈 东 高 


低 、 中 部 交错 分 布 的 状态 ,暴露 风险 人 口 密度 最 大 地 


级 市 共 16 个 ;地 下 水 储量 下 降 人 口 累积 百分比 呈现 明显 上 升 的 趋势 ,在 2016 年 达到 最 大 值 。 研 究 
结果 以 期 为 黄河 流域 地 下 水 资源 可 持续 利用 提供 科学 参考 。 
关 键 词 : GRACE; 地 下 水 储量 变化 ; 时 空 特征 ; 人 口 暴露 风险 ， 黄 河流 域 
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随 着 人 类 活动 加 强 , 水 资源 消耗 进一步 加 剧 ， 
人 水 矛盾 日 益 凸 显 "。 地 下 水 作为 水 资源 的 重要 组 
成 部 分 ,由 于 储量 稳定 ,水 质 好 ,可 被 各 行 各 业 直接 
使 用 ,是 人 类 活动 的 必需 品 汪 ,对 人 类 生存 发 展 与 
资源 可 持续 至 关 重 要 “。 然 而 ,部 分 地 区 由 于 不 合 
理 的 利用 ,导致 水 资源 问题 日 益 加 剧 ,引发 了 地 面 
沉降 ,海水 倒灌 、 地 下 水 资源 量 枯 竭 等 环境 问题 ,加 
剧 了 人 水 矛盾 ,尤其 是 在 干旱 半 干 旱地 区 更 为 突 
出 ”。2019 年 9 月 习近平 总 书记 在 河南 郑州 主持 
召开 黄河 流域 生态 保护 和 高 质量 发 展 座谈 会 并 发 
表 重 要 讲话 ,明确 指出 黄河 流域 在 我 国 经 济 社会 发 
展 和 生态 安全 方面 具有 十 分 重要 的 地 位 ,深刻 阐明 
黄河 流域 生态 保护 和 高 质量 发 展 的 重大 意义 ,作出 


收 稿 日 期 : 2022-03-04; ”修订 日 期 : 2022-05-19 


了 加 强 黄河 治理 保护 、 推 动 黄河 流域 高 质量 发 展 的 
重大 部 署 *。 在 此 背景 下 ,人 研究 黄河 流域 长 时 间 序 
列 的 地 下 水 储量 变化 ,对 该 区 域 的 水 资源 可 持续 利 
用 具有 重要 的 科学 意义 。 

传统 的 地 下 水 储量 监测 方法 耗 时 耗 力 , 且 在 监 
测 点 以 外 的 区 域 形成 了 监测 盲区 ,无 法 及 时 了 解 地 
下 水 储量 的 变化 情况 ”" ,而 重力 反 演 与 气候 实验 
(Gravity recovery and climate experiment, GRACE) E 
星 数 据 为 监测 大 尺度 区 域 地 下 水 储量 变化 提供 了 
有 效 途 径 ”。 国 内 外 众多 学 者 使 用 该 数据 开展 了 大 
量 研究 , 旦 主要 集中 在 地 下 水 储量 变化 监测 ,也 
有 部 分 学 者 开展 了 技术 方法 改进 .精度 检测 和 空间 
分 辩 率 降 尺 度 ”等 方面 的 研究 。 已 有 研究 在 量化 
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地 下 水 储量 的 动态 变化 方面 取得 丰富 的 成 果 , 对 本 
研究 具有 重要 的 借鉴 作用 ,但 鲜 有 文献 研究 地 下 水 
储量 变化 与 人 口 暴露 风 险 的 关系 ,而 这 方面 又 是 科 
学 制定 水 资源 政策 的 重要 依据 。 

黄河 流域 以 约 占 全 国 2% 的 径流 量 ,承载 了 全 
9% 的 人 口 ,水 资源 是 其 生产 .生活 的 制约 性 要 
素 , 故 本 文 以 黄河 流域 为 研究 区 ,基于 全 球 陆 面 数 
据 同化 系统 (Global land data assimilation system, GL- 
DAS) FI GRACE 遥感 数据 ,测算 了 黄河 流域 长 时 间 
序列 的 地 下 水 储量 变化 ,结合 LandScan 遥感 数据 ， 
探讨 地 下 水 储量 变化 与 人 口 暴 露 的 关系 , 主要 研究 
内 容 包括 :分 析 黄 河流 域 地 下 水 储量 变化 时 空 变化 
特征 、 识 别 地 下 水 储量 持续 下 降 的 区 域 , 并 分 析 该 
区 域 人 口 暴露 风险 情况 ,以 期 为 黄河 流域 地 下 水 资 
源 可 持续 利用 提供 科学 参考 。 


1 研究 区 概况 


黄河 发 源 于 青藏 高 原 ,流域 面积 达 75.24x10” 
km , 流 经 青海 省 、 四 川 省 .甘肃 省 .宁夏 回族 自治 
区 内 蒙古 自治 区 .陕西 省 .山西 省 .河南 省 .山东 省 
共 9 个 省 区 '( 图 1)。 黄 河 以 约 占 全 国 2% 的 径流 
量 ,承载 了 全 国 9% 的 人 口 ,另外 ,本 流域 和 流域 外 
引 黄 灌区 占 全 国 15% 的 耕地 面积 ,因此 ,黄河 流域 
在 中 国 经 济 社会 发 展 和 生态 安全 方面 具有 十 分 重 
要 的 地 位 。 由 于 其 主要 产 流 区 处 于 西北 干旱 半 干 


注 :该 图 基于 国家 测绘 地 理 信 息 局 标准 地 图 服务 网 站 下 载 的 审 
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旱 区 ,多 年 平均 降水 量 约 450 mm, 水 资源 量 少 且 时 空 
分 布 不 均 , 随 着 经 济 社会 的 不 断 发 展 , 该 区 域 取 水 量 
一 直 居 高 不 下 ,流域 水 资源 供需 矛盾 日 益 突出 中 。 
此 ,开展 黄河 流域 地 下 水 储量 变化 与 人 口 暴露 研 
究 ,对 于 黄河 水 资源 配置 战略 调整 ,实现 流域 生态 
保护 和 高 质量 发 展 具有 重要 意义 。 


2 数据 与 方法 


2.1 数据 来 源 与 处 理 

2.1.1 陆地 水 储量 数据 GRACE 卫星 上 搭载 了 精 
密 的 测 距 仪器 ,可 以 通过 双星 之 间 的 距离 变化 反 演 
地 球 表 层 局 部 重力 变化 情况 ,再 根据 重力 变化 推导 
出 监测 区 域内 陆地 总 水 储量 的 变化 "。 已 有 研究 
表明 ,基于 GRACE 卫星 数据 反 演 地 下 水 储量 与 监 
测 井 数 据 较 一 致 ,适用 于 中 国 区 域 的 地 下 水 储量 人 研 
究 。GRACE 数据 的 处 理 ,分 发 与 管理 由 Jet Propul- 
sion Lab( 美 国 喷气 动力 实验 室 ),GeoForschungsZen- 
trum (德国 地 学 研究 中 心 ) ,the University of Texas 
Center for Space Research( 美 国 德 克 萨 斯 大 学 空间 研 
究 中 心 )3 家 集体 负责 ,其 每 月 提供 地 球 重 力 场 模 
型 ,并 在 验证 后 分 发 给 公众 。GRACE 数据 产品 包括 
Level-0 ,Level-1A 、 Level-1B 和 Level-2 ,但 不 免费 对 
外 公布 ,为 了 方便 研究 质量 异常 (如 水 层 ) 的 用 户 ， 
一 些 团队 基于 Level-2 数 据 经 地 球 物理 校正 后 生成 
3 级 数据 产品 ,以 .nc 的 格式 存储 ,获取 方便 且 可 免费 


图 号 为 GS(2020)4619 号 的 标准 地 图 制作 , 底 图 边界 无 修改 。 下 同 。 


图 1 黄河 流域 概况 


Fig. 1 Overview of the Yellow River Basin 
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下 载 ,数据 链接 为 https://disc.gsfc.nasa.gov/datasets/ 
GLDAS_CLSM025_DA1_D_2.2/summary?keywords= 
CLSM%20eroundwater . 
2.1.2 土壤 含水 量 数据 土壤 含水 量 数据 来 源 于 戈 
达 德 地 球 科学 数据 与 信息 服务 中 心 (https://disc.gsfc. 
nasa.gov/datasets/GLDAS_NOAHO25_M 2.1/summary? 
keywords=GLDAS% 20CSRM) 。 该 研究 选取 了 GL- 
DAS-2.1 版 本 的 Noah 模型 输出 土壤 含水 量 数据 ,时 
间 序 列 为 2003 一 2016 年 的 逐 月 土壤 含水 量 数据 。 
数据 共 包含 4 层 土壤 水 分 数据 ,分别 为 地 下 0~10 
cm, Hb F 10~40 cm、 地 下 40~100 cm 和 地 下 100~200 
cm 共 4 层 ,土壤 水 分 数据 求 和 得 到 地 下 2 m 土壤 水 
分 。 将 2003 一 2016 年 每 月 土壤 含水 量 减 去 时 段 多 
年 平均 值 ,得 到 土壤 含水 量 距 平 值 , 便 可 得 到 逐 月 
土壤 含水 量变 化 。 
2.1.3 其 他 数据 ”人口 数据 来 源 于 LandScan 全 球 人 
口 格 网 数据 (https://landscan.ornl.gov) ,空间 分 辨 率 
为 1km, 已 有 大 量 学 者 使 用 该 数据 在 中 国 区 域 进 行 
精度 验证 与 科学 研究 “”。 植 被 覆盖 度 分 2 步 获 
得 :第 一 步 是 采用 最 大 值 合成 法 合成 年 NDVI 数 据 ， 
数据 获取 自 网 站 https:Wlandweb.nascom.nasa.gov/ ,时 
间 分 辩 率 为 16 d, 空 间 分 辨 率 为 230m。 第 二 步 利 
用 公式 求 得 植被 覆盖 度 。 
2.1.4 数据 预 处 理 首先 ,将 获取 所 有 原始 影像 数 
据 转 换 为 UTM 投影 ,地 理 坐 标 系 为 WGS84 坐标 ;其 
次 ,基于 黄河 流域 矢量 进行 掩 膜 处 理 ; 最 后 ,采用 重 
采样 工具 ,把 研究 所 需 的 所 有 栅 格 数据 重 采 样 至 最 
小 栅 格 分 辨 率 上 ,其 空间 分 辩 率 为 1 km 的 栅 格 数 
据 ,以 保持 栅 格 数据 像 元 大 小 的 一 致 性 ,时 间 分 辩 
率 以 月 为 尺度 进行 匹配 。 
2.2 研究 方法 
2.2.1 地 下 水 储量 计算 方法 ”根据 水 量 平衡 原理 计 
算 2003—2016 年 黄河 流域 地 下 水 储量 年 均 变 化 
量 。 陆 地 水 储量 的 变化 是 由 土壤 水 .地表 水 、 雪 / 冰 
川 融 水 以 及 地 下 水 储量 的 变化 引起 。 而 在 长 时 间 
尺度 上 ,土壤 含水 量 和 地 下 水 储量 变化 是 陆地 水 储 
量变 化 的 主要 原因 , 雪 水 当量 及 地 表 水 莫 水 量变 化 
可 以 忽略 不 计 "。 从 GRACE 反 演 的 陆地 水 储量 变 
化 中 减 去 土壤 含水 量变 化 得 到 地 下 水 储量 变化 。 
计算 公式 如 下 : 

AGWS, = ATWS, -ASMS， (1) 
式 中 :i 为 年 份 ;AGWS; 为 第 i 年 地 下 水 储量 年 均 变 化 


Œ (mm) ;ATWS :为 第 ;年 陆地 水 储量 年 均 变 化 量 
(mm); ASMS 为 第 ;年 土壤 水 年 均 变 化 量 (mm ) 。 
2.2.2 Theil-Sen Median #4 #4 4 Theil-Sen Medi- 
an 趋势 分 析 是 一 种 稳健 的 非 参 数 统计 趋势 计算 方 
法 ,可 用 来 反映 黄河 流域 地 下 水 储量 长 时 间 序 列 
演化 趋势 ,计算 公式 如 下 : 


{WW, es 
S, = Median( 1 } 2003 < i< j< 2016 (2) 


j-i 


SUES WASE j AE BS ee y i e F KAA E SK Te 
(mm); WJS i AE AO Se Te i Re FKE E SEROK 
(mm);Sv 反 映 黄河 流域 地 下 水 储量 等 效 水 高 变化 
的 单调 趋势 ,为 计算 zz-D/2 个 数据 组 合 的 斜率 中 
位 数 。 当 Sw>0 时 ,黄河 流域 地 下 水 储量 呈 增 长 趋 
势 ; 当 Svw<0 时 ,黄河 流域 地 下 水 储量 则 呈 减 少 趋势 。 
2.2.3 Mann-Kendall 检验 X H Mann- Kendall 方 
法 检验 黄河 流域 地 下 水 储量 时 间 趋 势 的 显著 性 , 设 
定时 间 序 列 {W} ,i=2003 2004, ++, 2016, HATE! 
如 下 : 


S-1 s>0 

s(S) 
定义 Z 统 计量 为 : Z=10， S=0 

S+1 S<0 (3) 

s(S) 


其 中 ,$= sen,-W) 


j=1i=j+1 


1, W,-W,>0 
senW,-W)=10, W-W.=0 (4) 
-1, W,-W,<0 

-H2 
ygy= nn en *>) (5) 


式 中 :5 为 变化 性 的 统计 变量 ;s(5) 为 5S 的 方差 ;n 为 时 
间 序 列 的 长 度 ;sgn( 田 为 符号 函数 。 当 IZl>pw, 表 
明 在 a 水 平 上 存在 显著 变化 。 该 研究 在 a=0.05 置信 
水 平 上 判断 黄河 流域 地 下 水 储量 时 间 序 列 上 变化 
趋势 的 显著 性 。 对 检验 结果 进行 分 级 ,按照 显著 性 
可 分 为 极 显 著 变 化 (1Z1>2.58) .显著 变化 (1.96<lZIs 
2.58) . 弱 显 著 变化 (1.65<lZ21< 1.96) 和 无 显著 变化 (0< 
1ZI<1.65)4 个 等 级 。 

2.2.4 地 下 水 储量 变化 分 级 方法 “参考 He GR 
用 标准 差分 级 法 ,根据 历年 地 下 水 储量 变化 均值 和 
标准 差 确定 每 段 分 割 的 端点 ,将 2003 一 2016 年 地 下 
水 储量 等 效 水 高 分 为 5 个 级 别 ,分 别 为 增加 、 稳 定 、 
减少 、 较 剧烈 减少 .剧烈 减少 ( 表 1)。 
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表 1 地 下 水 储量 变化 分 级 标准 


Tab.1 Grading standard for changes in groundwater storage 


地 下 水 储量 等 级 剧烈 减少 较 剧烈 减少 减少 稳定 增加 
分 级 标准 <x-1.5s XT—1.5s~x— 0.5s ¥-0.5s~¥+0.58 ¥+0.5s~¥+ 1.58 >¥+1.5s 
年 均 变 化 率 /mm'.an” <-58.38 —58.38~-34.54 -34.54~-10.70 -10.70~13.14 > 13.14 


TE: 元 为 中 国 2003 一 2016 年 地 下 水 储量 年 均 变 化 率 均 值 ;s 为 其 标准 差 。 


3 结果 与 分 析 


3.1 地 下 水 储量 空间 变化 特征 

2003 一 2016 年 黄河 流域 地 下 水 储量 年 均等 效 
水 高 在 空间 上 呈现 显著 的 东西 差异 ,东部 减少 西部 
增加 ,多 年 平均 等 效 水 高 为 -22.62 mm- a", FH PG IH 
东 地 下 水 储量 的 下 降 程度 不 断 加 剧 , 且 下 降 区 域 呈 
现 由 下 游 向 中 上 游 扩 散 的 态势 (图 2)。 其 中 ,下 降 
区 域 主要 集中 在 中 东部 地 区 ,变化 等 级 以 较 剧 烈 减 
少 和 剧烈 减少 为 主 。 研 究 初期 ,剧烈 减少 的 区 域 包 
含 西安 市 .洛阳 市 三门峡 市 .郑州 市 等 城市 ;研究 
末期 ,已 包含 黄河 中 下 游 的 大 部 分 城市 ,甚至 上 游 


游 地 区 经 济 增 长 快 . 人 口 集聚 能 力 强 ,使 得 农业 T 
业 和 生活 用 水 不 断 增加 ,而 地 下 水 补充 速率 远 远 低 
于 消耗 速率 ,故地 下 水 储量 不 断 减少 。 而 西南 部 的 
地 下 水 储量 仍 以 增加 为 主 ,这 是 由 于 该 区 域 的 地 形 
以 青藏 高 原 和 黄土 高 原 为 主 ,平均 海拔 2000 m 以 
上 ,人口 稀少 且 工 农业 欠 发 达 , 地 下 水 资源 量 消耗 
总 量 少 ,加 之 全 球 气候 变化 引起 局 部 地 区 冰川 融 水 
增加 ,使 得 黄河 流域 上 游 地 下 水 补给 速率 大 于 消耗 
速率 ,因此 地 下 水 储量 变化 呈 增 加 趋势 。 

3.2 地 下 水 储量 时 序 变化 特征 

3.2.1 RRR HE 通过 分 区 统计 黄河 流域 
2003 一 2016 年 地 下 水 储量 变化 情况 ,由 图 3 A 
下 水 储量 变化 整体 呈 下 降 趋势 。 该 区 域 地 下 水 储 


的 玉树 藏族 自治 州 . 果 洛 藏族 自治 州 .海南 藏族 自 
治 州 等 也 成 为 剧烈 减少 区 域 。 这 是 由 于 黄河 中 下 


(g) 2009 年 


(k) 2013 年 


量 等 效 水 高 最 高 值 为 2004 年 的 10.32 mm, 最 低 值 为 
2016 年 的 -75.56 mm, 2003—2016 年 地 下 水 储量 年 


À () 2014 年 


>z 


0 200 km 0 200 km 
(0) 2013 一 2016 年 均值 
i 地 下 水 储量 等 效 水 高 分 级 
m 剧烈 减少 
il JD 

= e 图 例 
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图 2 2003 一 2016 年 黄河 流域 地 下 水 储量 空间 分 布 
Fig. 2 Spatial distribution of groundwater storage in the Yellow River Basin from 2003 to 2016 
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图 3 2003 一 2016 年 黄河 流域 地 下 水 储量 与 植被 覆盖 度 
变化 关系 


Fig. 3 Relationship between groundwater storage and 


fractional vegetation cover in the Yellow River Basin from 
2003 to 2016 


均 降 幅 为 5.93 mma, FLA, 2003—2007 年 黄河 流 
域 地 下 水 储量 等 效 水 高 变化 幅度 较 小 ,至 2008 年 
时 ,黄河 流域 地 下 水 储量 主要 变化 特征 由 增加 变 为 
减少 ,而 2009 一 2016 年 则 呈现 出 持续 下 降 的 趋势 ， 
特别 是 2015 一 2016 年 降幅 最 大 ,下 降 24.41 mm, 最 
大 降幅 为 47.74%。 

统计 2003 一 2016 年 黄河 流域 地 下 水 储量 历年 
春季 (3 一 5 月 )、 夏 季 (6 一 8 月 ) .秋季 (9 一 11 月 ) 和 冬 
季 (12 一 次 年 2 月 ) 的 均值 ,发 现 呈 现 明显 季节 变化 
特点 (图 4)。 春 冬季 下 降幅 度 大 ,夏秋 季 下 降幅 度 


减 小 。 这 是 因为 黄河 流域 受 季 风气 候 影响 ,夏秋 季 
是 雨季 ,地 表 径 流 补给 地 下 水 ,而 冬 春季 为 枯 季 ,地 
下 水 补给 地 表 水 。 地 下 水 储量 下 降幅 度 春 季 最 大 、 


冬季 次 之 、 秋 季 再 次 之 .夏季 最 小 。 其 中 ,春季 地 下 
水 储量 等 效 水 高 均值 为 -58.92 mm, 冬 季 地 下 水 储 
量 等 效 水 高 均值 为 -54.12 mm, 秋季 地 下 水 储量 等 
效 水 高 均值 为 -47.81 mm ,夏季 地 下 水 储量 等 效 水 
高 均值 为 -38.02 mm, 
3.2.2 时 序 趋势 将 征 ”为 了 进一步 刻画 黄河 流域 地 
下 水 储量 时 间 序 列 的 线性 变化 趋势 ,利用 Theil-Sen 
Median 和 Mann-Kendall 方 法 将 S$S 值 与 ZI 值 重 分 类 后 
再 进行 同 像 元 至 加 ,得 到 黄河 流域 地 下 水 储量 线性 
显著 性 变化 趋势 (图 5)。 黄 河流 域 地 下 水 储量 线性 
显著 性 变化 趋势 共 分 为 7 个 类 别 ,分 别 是 弱 显 著 减 
少 、 显 著 减 少 极 显著 减少 、 弱 显著 增加 、 显 著 增 加 、 
极 显 著 增 加 和 无 显著 变化 ,各 类 面积 见 表 2。 结 合 
图 5 与 表 2 可 知 , 黄 河流 域 地 下 水 储量 时 序 线性 显 
著 性 变化 特征 以 减少 为 主 。 其 中 极 显著 减少 的 区 
域 主要 分 布 在 中 下 游 地 区 ,包括 陕西 省 中 北部 和 山 
西 省 大 部 分 地 区 ,河南 省 和 山东 省 的 黄河 沿岸 地 区 
等 ,总 面积 为 50.62x10:km2, 占 黄河 流域 总 面积 的 
63.60% ; 显著 减少 的 区 域 分 布 在 青海 省 海北 藏族 自 
治 州 和 西宁 市 ,总 面积 为 2.36x10' km2, 占 黄河 流域 
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at 
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图 例 一 省 界 市 界 
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图 4 黄河 流域 地 下 水 储量 四 季 变 化 


Fig. 4 Seasonal variation of groundwater storage in the Yellow River Basin 
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图 5 黄河 流域 地 下 水 储量 线 履 


E 显 著 性 变化 趋势 


Fig. 5 Linear significant change trend of groundwater storage in the Yellow River Basin 


表 2 黄河 流域 地 下 水 储量 面积 变化 


Tab. 2 Variation of groundwater storage area in the Yellow River Basin 


变化 类 型 弱 显 著 减 少 。 显著 减少 极 显 著 减少 弱 显 著 增 加 显著 增加 极 显著 增加 无 显著 变化 
面积 /10:km?” 15.65 2.36 50.62 3.13 0.00043 6.13 1.69 
百分比 /% 19.66 2.96 63.60 3.94 0.00050 7.70 2.13 


总 面积 的 2.96%。 综 上 ,黄河 流域 地 下 水 储量 主要 
时 序 变 化 特征 为 持续 减少 ,这 是 由 人 文 与 自然 因素 
共同 作用 形成 ,一 方面 工农 业 生产 和 生活 用 水 总 量 
在 不 断 增 加 ,人 水 矛盾 愈 发 突出 ; 另 一 方面 ,生态 用 


地 下 水 暴露 风险 区 。 将 暴露 风险 区 人 口 密 度 分 为 6 
个 级 别 忠 ,分别 为 无 分 布 .1~20 人 :km?、21~50 人 … 
km- 51~100 人 .km 、 101~200 人 km” 了 和 >200 人， 
km”。 以 2016 年 为 例 ,统计 了 人 研究 区 内 市 级 行政 单 


水 在 逐年 加 大 ,如 研究 发 现 植 被 变化 与 地 下 水 储量 
旺 显著 的 负 相 关 关 系 。2003 一 2016 年 黄河 流域 植 
被 覆盖 度 时 间 序 列 上 的 变化 趋势 呈现 显著 的 上 升 


元 地 下 水 暴露 风险 下 人 口 密度 (图 6)。 黄 河流 域 地 
下 水 暴露 风险 状态 下 的 地 级 市 人 口 密度 在 空间 上 
呈 东 高 西 低 .中 部 交错 分 布 的 状态 。 其 中 ,暴露 风 


趋势 ,与 地 下 水 储量 的 年 际 变化 趋势 相反 。2003 年 
植被 覆盖 度 为 0.30,2016 年 增长 至 0.35 ,时 稳步 增长 
的 状态 ,地 表 植 被 覆盖 情况 日 益 改 善 (图 3)。 使 用 
SPSS 对 两 者 进行 相关 性 分 析 , 发 现在 P<0.01 显 著 性 
水 平 上 ,皮尔 逊 相关 系数 为 -0.97, 呈 显著 的 负 相 关 
关系 。 黄 河流 域 气 候 类 型 以 干旱 半 干 旱 气 候 为 主 ， 
地 下 水 储量 与 植被 生态 环境 的 关系 相 较 于 其 他 气 
候 区 更 为 密切 ,植被 覆盖 度 的 增加 会 消耗 大 量 的 水 
资源 , 且 以 地 下 水 为 主 , 说 明 植 被 覆盖 度 变化 对 当 


险 区 人 口 密度 >200 人 .km 的 地 级 市 共 16 个 ,包括 
兰州 市 .银川 市 西安 市 渭南 市 . 晋 中 市 .运城 市 、 
临汾 市 济源 市 .安阳 市 , 浴 阳 市 .焦作 市 .新 乡 市 、 
郑州 市 .菏泽 市 .泰安 市 和 济南 市 ;暴露 风险 区 人 口 
密度 >100 人 :km 上 且 <200 人 :ko 的 地 级 市 共 5 个 ， 
分 别 为 临夏 回族 自治 州 `. 乌 兰 察 布 .晋城 市 吕梁 市 
和 洛阳 市 ;暴露 风险 区 人 口 密度 >50 A km” H< 100 
人 :km ”的 地 级 市 共 4 个 ,分 别 为 白银 市 平凉 市 . 包 


地 地 下 水 水 资源 储量 产生 深刻 影响 。 

3.3 地 下 水 储量 下 降 区 人 口 暴露 风险 状态 

3.3.1 暴露 风险 状态 下 人 口 分 布 为 精确 识别 黄河 
流域 地 下 水 暴露 风险 区 ,将 研究 区 地 下 水 储量 线性 
变化 趋势 中 显著 减少 与 极 显著 减少 的 区 域 定 义 为 


头 市 和 三 门 峡 市 ;暴露 风险 区 人 口 密度 为 0~50 人， 
km” 的 地 级 市 交错 分 布 。 

3.3.2 地 下 水 储量 下 降 人 口 占 比 ”为 便于 分 析 , 将 
2004 一 2016 年 每 隔 4a 统 计 地 下 水 储量 等 效 水 高 变 
化 幅度 分 为 17 个 不 等 间距 的 区 间 ,并 计算 超过 各 级 
限 值 的 人 口 累 积 百分比 (图 7)。 结 果 表 明 超 过 各 级 
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图 6 黄河 流域 地 下 水 暴露 风险 下 人 


Fig.6 Spatial distribution of population density under groundwater exposure risk in the Yellow River Basin 
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7 黄河 流域 地 下 水 储量 的 下 降水 平 超过 各 级 限 值 的 人 
口 累积 百分比 
Fig. 7 Percentage of population whose groundwater storage 
decline level exceeds the limits at all levels in the Yellow 


River Basin 


限 值 的 暴露 人 口 累积 百分比 曲线 在 2004 一 2016 年 
向 右 移动 , 呈 明显 上 升 的 趋势 。 其 中 ,等 效 水 高 <0 
mm 的 区 域 人 口 在 2004 年 有 1517.62x10' 人 , 占 黄河 
流域 总 人 口 的 25.80% , 2008 4E _E Ft Æ 3705.02 10* 
人 , 占 比 上 升 至 61.54% ,2012 年 占 比 达到 82.55% , 
2016 年 占 比 达到 最 高 为 87.04%, 地 下 水 储量 下 降 区 
总 人 口上 升 至 3908.94x10' 人 ;等 效 水 高 <-15 mm 
的 区 域 人 口 在 2004 年 有 445.96x10'* 人 , 仅 占 区 域 
总 人 口 的 7.58%,2008 年 陡 升 至 31.26%,2012 年 更 
是 达到 73.08% ,2016 年 上 升 至 最 高 79.07% ;等 效 
水 高 <-35 mm 的 区 域 人 口 在 2004 年 有 16.46x10 
人 , 仅 占 流域 总 人 口 的 0.28% ,2008 年 上 升 至 
2.60% ,2012 年 骤然 上 升 至 63.12% ,2016 年 高 达到 
66.35% ;等 效 水 高 <-55 mm 的 区 域 人 口 ,2012 占 区 


口 密度 空间 分 布 


域 总 人 口 的 48.26% ,2016 年 占 比 达 60.72% ;等 效 水 
高 <-65 mm 的 区 域 人 口 在 2012 年 占 区 域 总 人 口 的 
34.17% ,2016 年 上 升 为 57.4496。 


4 讨论 


地 下 水 储量 急剧 下 降 不 仅 威胁 到 城市 居民 的 
用 水 安全 ,还 会 引起 地 面 沉降 等 问题 ,成 为 威胁 区 
域 城市 发 展 及 居民 安全 的 重大 灾害 隐患 的 。 目 前 ， 
干旱 半 干 旱 区 地 下 水 储量 长 时 间 序 列 宏观 尺度 的 
人 研究 还 较 少 ,这 是 由 于 中 国 地 质 结构 复杂 ,含水 层 
较 破 碎 3, 且 地 下 水 储量 人 工 监测 成 本 高 .难度 
大 。 该 研究 融合 GRACE、GLDAS 和 LandScan A O 
数据 等 多 种 遥感 数据 ,并 参考 前 人 研究 思路 ,根据 
流域 水 循环 与 水 量 平衡 原理 ,测算 黄河 流域 2003 一 
2016 年 地 下 水 储量 ,识别 其 时 空 变化 趋势 与 特征 ， 
分 析 该 区 域 地 下 水 储量 在 不 同 降级 下 的 人 口 暴露 
特征 。 

研究 结果 表明 黄河 流域 地 下 水 储量 逐年 下 降 ， 
人 口 暴露 占 比 持续 上 升 。 对 比 已 有 人 研究 ,如 冯 伟 
等 9 得 出 华北 平原 2002 一 2014 年 地 下 水 等 效 水 高 
质量 亏损 为 (56+6)mm*a'!; Yi 等 ?发现 2003 一 2014 
年 黄河 -淮河 -辽河 -海河 流域 的 地 下 水 质量 损失 也 
比较 明显 ,损失 的 质量 为 (10.2+0.9)Gt.a'; 李 婉 秋 
等 5 研究 表明 关中 地 区 2003 一 2014 年 地 下 水 储量 
存在 长 期 亏损 ,亏损 速率 为 3.70 mma'。 综 上 ,本 
研究 结论 与 已 有 研究 基本 一 致 。 基 于 此 ,在 黄河 流 
域 生态 保护 与 高 质量 发 展 的 背景 下 ,相关 部 门 应 将 
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确 地 探究 干旱 半 干 旱 区 地 下 水 时 空 动态 变化 过 程 


应 积极 发 展 节 水 农业 .重点 区 域 采 用 淡水 回 灌 补充 
地 下 水 等 手段 ,保证 黄河 流域 地 下 水 储量 安全 ™。 


与 影响 机 制 。 第 二 ,黄河 流域 地 下 水 储量 变化 主要 
集中 在 中 下 游 地 区 ,是 本 文 的 研究 重点 ,而 上 游 高 


该 研究 不 但 揭示 了 地 下 水 储量 变化 的 时 空 特征 ,还 
探究 了 流域 不 同 地 区 人 口 暴露 分 析 的 空间 异 质 性 ， 
研究 结果 有 助 于 厘清 黄河 流域 人 水 矛盾 的 重点 区 
域 ,对 实现 流域 可 持续 发 展 具有 重要 意义 。 

然而 ,本 研究 可 能 存在 以 下 不 足 。(1) 数据 方 
面 。 首 先 ,考虑 数据 可 获得 性 ,本文 只 选取 一 代 
GRACE 数据 ,时 间 序 列 为 2003 一 2016 年 ,这 就 限制 
研究 的 时 效 性 。 同 时 ,本 研究 所 用 的 GRACE 数据 ， 
虽 经 过 了 数据 校正 和 地 理 信息 技术 处 理 ,但 是 受制 
于 时 空 分 辩 率 限制 ,仍然 存在 一 定 的 不 确定 性 。 此 
外 ,偶然 因素 也 影响 GRACE 部 分 月 数据 精度 , 李 杰 
等 中 对 GRACE 时间 序 列 进行 月 趋势 项 和 平均 值 ， 
孙 倩 等 ”采用 三 点 滑动 平均 法 消除 了 季节 内 的 随 
机 波动 。 此 外 ,对 于 土壤 水 分 数据 的 处 理 也 有 无 法 
避免 的 不 确定 性 ,因为 许多 地 区 地 表 水 文 过 程 较为 
复杂 且 地 表 水 数据 缺乏 ,许多 研究 计算 所 得 的 地 表 
水 储量 变化 与 实际 值 往往 相差 较 大 。 为 减 小 计算 
地 表 水 数据 过 程 中 的 不 确定 性 ,土壤 含水 量 数据 由 
GLDAS-NOAH 模 型 输出 ,而 NOAH 模 型 中 只 考虑 了 
4 层 土壤 含水 量 , 仅 包含 地 表 下 2 m 的 土壤 含水 量 。 
为 此 ,有 部 分 学 者 计算 了 GLDAS 的 CLM MOSAIC, 
NOAH Fil VIC 这 4 个 模型 输出 的 土壤 含水 量 数据 ， 
分 别 对 应 3.43 m、3.50 m,2.00 m Ail 1.90 m 土壤 舍 水 
量 ,然后 取 其 平均 值 。 尽 管 本 人 研究 在 充分 参考 前 人 
的 理论 方法 ,最 终 确 定 选 取 GLDAS-NOAH 模 型 输出 
的 土壤 含水 量 数据 ,包含 地 表 下 2m 的 土壤 水 含量 ， 
但 仍 存在 一 定 不 确定 性 。(2) 对 黄河 上 游 地 下 水 储 
量 计算 有 待 进一步 细 化 ,尽管 本 研究 中 借鉴 前 人 思 
路 ,认为 土壤 含水 量 和 地 下 水 储量 变化 是 陆地 水 储 
量变 化 的 主要 原因 , 雪 水 当量 及 地 表 水 蓄 水 量变 化 
可 以 忽略 不 计 史 ,这 对 于 黄河 中 下 游 是 有 效 的 ,但 
若 重点 关注 黄河 上 游 地 区 地 下 水 储量 的 变化 ,充分 


海拔 地 区 的 地 下 水 储量 变化 小 ,人 类 经 济 活 动 对 其 
影响 较 小 ,并 考虑 到 GRACE 数据 的 空间 分 辩 率 问 
题 ,因此 在 黄河 流域 地 下 水 储量 变化 计算 时 只 考虑 
了 土壤 水 分 的 影响 ,下 一 步 工作 中 应 补充 雪 水 和 植 
物 水 变化 对 地 下 水 储量 变化 的 影响 。 第 三 ,数据 方 
面 ,结合 便携 式 绝对 重力 仪 和 地 下 水 监测 井 数据 对 
地 下 水 储量 的 空间 分 布 及 变化 趋势 进行 反 演 、 验 证 
和 分 析 , 以 便 进一步 研究 人 水 矛盾 突出 区 域 的 地 下 
水 资源 演变 规律 ,为 黄河 流域 未 来 经 济 发 展 与 资源 
可 持续 利用 提供 科学 参考 。 


5 结论 


(1) 从 空间 特征 看 ,黄河 流域 地 下 水 储量 变化 
特征 在 空间 上 旦 明显 的 东西 差异 ,东部 减少 西部 增 
加 ,多 年 平均 等 效 水 高 为 -22.62 mm: a`, H P A 
地 下 水 储量 的 下 降 程 度 不 断 加 剧 , 且 下 降 区 域 呈 由 
下 游 向 中 上 游 扩散 的 态势 。 其 中 ,下 降 区 域 主要 集 
中 在 中 东部 地 区 ,变化 等 级 以 较 剧 烈 减少 和 剧烈 减 
少 为 主 。 研 究 初期 ,剧烈 减少 的 区 域 包含 西安 市 、 
洛阳 市 三门峡 市 .郑州 市 等 城市 ;研究 未 期 ,已 包 
含 黄河 中 下 游 的 大 部 分 城市 ,甚至 上 游 的 玉树 藏族 
自治 州 、 果 洛 藏族 自治 州 海南 藏族 自治 州 等 也 成 
为 剧烈 减少 区 域 。 表 明黄 河流 域 人 水 关系 矛盾 日 
趋 突 出 ,该 区 域 经 济 发 展 一 定 程 度 上 依赖 于 地 下 水 
资源 的 过 度 开 采 。 

(2) 从 时 间 特 征 看 ,一 方面 ,2003 一 2016 年 黄河 
流域 地 下 水 储量 时 间 变 化 特征 为 下 降 。 该 区 域 地 
下 水 等 效 水 高 最 高 值 为 2004 年 的 10.32 mm ,最低 值 
为 2016 年 的 -75.56 mm,2003 一 2016 年 地 下 水 储量 
年 均 降 幅 为 5.93 mma'。 其 中 ,2015 一 2016 年 降幅 
最 大 ,下 降 24.41 mm, 最 大 降幅 为 47.74%。 男 一 方 
面 , 地 下 水 储量 时 间 变 化 存在 季节 效应 ,地 下 水 储 


考虑 雪 水 .植物 根 区 水 和 冠 层 水 等 对 地 下 水 储量 的 
影响 是 非常 有 必要 的 。(3) 地 下 水 储量 变化 与 人 类 
经 济 活动 的 关系 有 待 进一步 挖掘。 

在 未 来 的 工作 中 ,以 下 3 个 方面 有 待 进一步 研 
究 。 第 一 ,基于 CRACE 数 据 , 结 合 地 下 水 实测 站 点 
数据 , 反 演 中 国 更 高 时 空 分 辨 率 地 下 水 储量 数据 ， 
并 填补 两 代 GRACE 卫星 数据 空白 期 数据 ,以 更 精 


量 下 降幅 度 春 季 最 大 、 冬 季 次 之 秋季 上 青 次 之 .夏季 
最 小 。 春 季 地 下 水 等 效 水 高 均值 为 -58.92 mm, & 
季 地 下 水 等 效 水 高 均值 为 -54.12 mm, ,秋季 地 下 水 
等 效 水 高 均值 为 -47.81 mm ,夏季 地 下 水 等 效 水 高 
均值 为 -38.02 mm. 

(3) 从 人 口 暴 露 风险 看 ,黄河 流域 地 下 水 暴露 
风险 状态 下 的 地 级 市 人 口 密度 在 空间 上 呈 东 高 西 
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低 .中 部 交错 分 布 的 状态 ,暴露 
大 地 级 市 共 


风险 区 人 口 密度 最 
16 个 ,包括 兰州 市 .银川 市 .西安 市 、 涓 


南 市 
过 各 级 限 值 的 人 口 累积 


Pi .运城 市 等 。 地 下 水 储量 下 降幅 度 超 
分 比 不 断 增加 ,在 2016 年 


达到 最 大 值 ,人 水 矛盾 十 分 突出 。 
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Groundwater storage and population exposure in the Yellow River Basin 
based on GRACE data 


DENG Chun", JIANG Xiaohui, SUN Weifeng’ 


(1. College of Urban and Environmental Sciences, Northwest University/Shaanxi Key Laboratory of Earth Surface System and 
Environmental Carrying Capacity, Xi an 710127, Shaanxi, China; 2. Department of Economics and Management, Yuncheng 


University, Yuncheng 044000, Shaanxi, China) 


Abstract: In recent years, the contradiction between humans and water resources in the Yellow River Basin, 
China has become increasingly prominent, particularly the excessive consumption of groundwater has become 
one of the main contradictions restricting the harmonious development of humans and nature in this region. 
Based on GRACE satellite data and GLDAS data, the change in groundwater storage of the Yellow River Basin 
was calculated from 2003 to 2016, then its temporal and spatial characteristics were explored, and population 
exposure in areas with significant groundwater decline was identified. The results show the following: (1) In 
terms of spatial distribution, the groundwater storage of the Yellow River Basin are more in the west and less in 
the east, the degree of decline from the west to east continues to intensify, and the decline area shows a spreading 
trend from the lower reaches to the middle and upper reaches. (2) The groundwater storage showed a downward 
trend in the Yellow River Basin from 2003 to 2016, with an average annual decrease of 5.93 mm-a '. It is a 
continuous decline from 2004 to 2016, and the biggest drop occurs in 2015—2016. In addition, there is a seasonal 
effect on the time change of groundwater storage, groundwater storage declines the most in spring, the second 
most in winter, the third most in autumn, and the least in summer. (3) Based on the risk of groundwater exposure, 
the characteristics of population density of prefecture-level cities are high in the east, low in the west, and have a 
staggered distribution in the middle, and there are 16 prefecture-level cities with the highest population density at 
risk of exposure. The cumulative percentage of the population with groundwater storage decline exhibited a clear 
upward trend, reaching a maximum in 2016. The results provide scientific reference for sustainable utilization of 
groundwater resources in the Yellow River Basin. 

Key words: GRACE; groundwater storage change; spatial-temporal characteristics; population exposure risk; 
Yellow River Basin 


